Electron localization/delocalization in incommensurate helical magnets 
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The electronic states in incommensurate (IC) helical magnets are studied theoretically from the 
viewpoint of the localization/delocalization. It is found that in the multi-band system with rela- 
tivistic spin-orbit interaction, the electronic wavefunctions show both the extended and localized 
natures along the helical axis depending on the orbital, helical wavenumber, and the direction of 
the plane on which spins rotate. Possible realization of this localization is discussed. 
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The helical magnets have been studied for a long term 
since its first discovery by Yoshimori Jj. Their ground 
states are determined by the (frustrated) exchange inter- 
actions and their Fourier transformation J{q). Various 
properties including the spin wave excitations are ana- 
lyzed theoretically for many materials. (See 0] for an 
early review.) Helical spin structure is now attracting 
recent revived interests from the viewpoints of both di- 
electric and transport properties. One example is the 
ferroelectricity induced by the helical magnetic order. 
Theoretically, the spin current associated with the non- 
collinear spin configuration is proposed to induce the 
electric polarization Q. Experimentally, it is now found 
that this mechanism is at work in fJMnOa ^ |^, |^ and 
in other materials 0, H, IS Q • Another new aspect is 
the anomalous transport properties associated with the 
onset of helical spin structure in metallic systems such 
as/3-Mn02 SrFeOa [H [3 and MnSi, (Fe,Co)Si 

|l4|. These developments urge the microscopic theory 
of electronic states to understand the physical properties 
associated with the helical spins. 

In the absence of the spin-orbit interaction (SOI), one 
can rotate spin frame so that the z-axis is parallel to the 
direction of the local spin. In this rotated frame, the 
spins are aligned ferromagnetically and the original spin 
structure is reflected in the magnitude and phase of the 
effective transfer integrals. This leads to the double ex- 
change interaction 1151 and various phenomena related 
to the spin chiralitv |l6l Il7l | , respectively. When we con- 
sider the state of single helical wavevector q, the relative 
angle between the neighboring spins does not break the 
original translational symmetry. Furthermore, there is 
no spin chirality, i.e., no fictitious magnetic field induced 
by the solid angle subtended by the spins. Therefore the 
Hamiltonian in the rotated spin frame preserves the pe- 
riodicity of the original lattice, and hence one can define 
the Bloch wavefunction. 

This situation is modified in an essential way when the 



SOI is taken into account. In this case, one cannot ro- 
tate the spin frame with the orbitals being intact, and 
the transfer integrals forming a matrix between ions are 
transformed in a nontrivial way. Therefore, in general, 
we expect the incommensurate (IC) modulation of the 
transfer integrals and even of the site energies in the ef- 
fective Hamiltonian in the rotated frame. 

Localization/delocalization of electronic state in an IC 
potential is an old issue Unlike in the case of com- 
mensurate periodic potentials, the eigen states are not 
the extended Bloch states in the case of IC potentials. 
Therefore the band structures would be unusual, i.e., 
highly fragmented, in those IC potentials. The central is- 
sue is whether electronic states are extended or localized 
in such kind of potentials, namely, Metal-Insulator Tran- 
sition (MIT). Aubry and Andre (A-A) 19] have shown 
that in a simple ID model MIT occurs simultaneously for 
all energies when the strength of the IC potential Vb is 
equal to the transfer integral t, i.e., if Vo is greater than 
t, the electronic states localize. We can also regard A-A 
model as a two dimensional tight binding model with IC 
magnetic flux. Actually, well-known Hofstadter butterfly 
is closely related to this model (2fli | . Using the trace map 
technique, Kohmoto et al. [i^ has exactly studied the 
scaling properties of the Fibonacci lattice system which 
can be regarded as the A-A model with IC modulation 
Qa/2n being the inverse Golden Mean. Similar prob- 
lems with IC transfer integral are also investigated by 
Kohmoto et al. Ej ]. 

In this paper, we investigate the localiza- 
tion/delocalization of electronic states in IC helical 
magnets. First we study a model of 5 c?-orbitals in 
cubic symmetry taking into account the SOI. We found 
that as SOI increases, the localization caused by IC 
starts from the specific t2g wavefunctions at around 
q ^ ir/a {a: lattice constant). In order to scrutinize 
this localization, we construct an effective single-band 
model for t2g bands. With this effective model, the 
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localization lengths are studied in more detail including 
its dependence on the angle ip between the spin rotation 
angle and the helical wavevector. 

We start with the following electronic model: 

H = Hu + Hso + Hd + Ht, 

3 

Hso = -A^ • S;j, 

3 

Hd - ^e„M^J(d^J, 

3 

Ht = tf^\d'^,){d'^+,J. (1) 

In the octahedral ligand field, the d orbitals are split into 
Cg- and f2g-orbitals H^. The f2g-orbitals, i.e., d^y,dy'^ 
and d^^ , have energies lower than e^-orbitals, i.e., d^^"^^, 

r, 2 2 

and c? by lOZJg, but the order is reversed as we take 

the hole picture in the foUowings, i.e., et^g ~ = lODq. 
The on-site SOI is considered, the matrix elements of 
which are calculated by L- S with L {S) being the orbital 
(spin) angular momentum. It is noted that L has no ma- 
trix elements within the Cg sector, while nonzero coupling 
occurs within t2g sector and between Cg and t2g sectors. 
Conside ring: the hopping between d-orbitals and oxygen 
orbitals |2J|, we derive the effective transfer integrals tfj 
between d"-orbitals at neighboring magnetic ions i and 
j. We took the values t^^ = = 0.1, i^^""''" = 0.3 
and = f^^-y^ = 0. In Hu of Eq. (1), the mag- 
netic moment at site j is described by the unit vector 
Cj = {cos (pjSinO j ,sin(j)jsm6j , cos 9j) and Sj denotes the 
electronic spin operator at site j. We assume the IC 
helical magnetic structure for Sj along 2;-axis, which is 
on the spin {z, a;)-plane, realized as a result of the frus- 
trated spin exchange interaction. We focus on the or- 
dered ground state properties, and hence the mean field 
treatment gives a good description of the system. We 
assume the ferromagnetic spin configuration perpendic- 
ular to the helical wavevector q, and hence k^, ky are 
good quantum numbers, i.e., the electronic wavefunc- 
tions are plane waves along x- and y-directions. We fix 
kx — ky = hereafter, and consider the one-dimensional 
(ID) model only along z-direction. Fig. 1 shows the 
calculated density of states as a function of the helical 
wavenumber q with the color specifying the localization 
length ^. We note here that the sample size is a prime 
number 199, and helical wavenumbers g's are taken to be 
proximate to the IC values. All the bands states from the 
Cg-orbitals are extended due to the weak SOI, and hence 
are omitted in Fig. ^ The green region is the extended 
states while the blue one is strongly localized within the 
scale of lattice constant. We took the values lODq — 3, 
U = 1.4, and A = 1.0. 

The density of states are understood as follows. The 
largest splitting between Cg and t2g occurs due to the lig- 
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FIG. 1: Density of states and color map of the localization 
length 4 on it for the d-orbital model Eq. (1). 



and field lODq in Eq. (1). Then the t2g bands are further 
split into bands of Fy and Fg origin, the latter of which 
is upper in energy since we take the hole picture. Then 
both the bands are further split by the spin exchange 
field U. 

By using the iterative method developed by MacKin- 
non [2^, we can calculate the Green's function = 
{1\{E — H)~^\N), which connects both ends of the long 
chain. is still a 10 x 10 matrix and the Lyapnov 

exponent, i.e., the inverse of the localization length ^ 

is obtained as | = — limTv^oo 2W^'^'^'^l^i^l^- "^^^ ^^^^ 
color at around q ^ n/a means the strong localization 
along the helical axis. When we change A, we still ob- 
serve the localization down to A ~ 0.2. Therefore we 
conclude that the localization starts in some part of the 
electronic spectrum at around q = n/a as one increases 
the SOI. The most remarkable point we can grasp from 
the above figure is the energy dependence of localiza- 
tion/delocalization nature, namely that there are both 
localized and extended states at different energies for the 
same q. This is in sharp contrast to the case of A-A model 
where all the states are either extended or localized de- 
pending only on the ratio Vo/t as mentioned above. 
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In order to study this localization in more depth, we 
now derive the effective model for a limiting case, i.e., 
lODq :s> 3A/2 ^ U ^ t. Even though this is not nesces- 
sarily a suitable limit for realisitc systems, it clarifies why 
^ depends on the orbitals. By taking into account the 
spin degree of freedom, there is six-fold degeneracy of 
the t2g energy levels. Because of the on-site SOI, this 
degeneracy is lifted and we have t wo g roups of spin-orbit 
coupled states, labeled r7 and Fg |23|. The two-fold de- 
generate states, Fy, and the four- fold degenerate one, Fg, 
are given by |3+) = dd^'') + |df ) + ))/V3, 13") = 

)-|-z|df ))/y3, and |1+) = (|df )+i|df ))/\/2, 
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= (|df ) - ))/V2, |2+> ^ - \df) - 

iMf ))/V6, |2-) = (2|d7) + |df ) - iMf ) )/V6, respec- 
tively, where the quantization axis of spin is taken to be 
the z axis. Henceforth, we assume that the spin-orbit 
coupling in our system is sufficiently large and focus only 
on the case where the two multiplets, i.e., Fy and Fg, do 
not hybridize with each other. 

Now, we construct the normalized state \gj) to mini- 
mize {gj \ — Ucj ■ Sj\gj) in the Hilbert space spanned by 
the states in F7 or Fg. The desired states whose spins 
are parallel to the unit vector ej are explicitly given for 
F7 and Fg by 



cos^|3- 



(2) 



the two spins Sn and S„+i. The phase a„_„+i is the 
vector potential generated by the noncoUinear spin con- 
figuration, but we can eliminate it by appropriate gauge 
transformation. Then we can conclude that we have no 
incommensurability in our ID F7 model. 




FIG. 2; Spin plane tilted by angle (p from xy-plane (left). 
The helical spins are rotating on the tilted plane placed peri- 
odically placed along the z-axis. Blue arrows represent spins, 
while a denotes the lattice spacing (right). 
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respectively. Here, subscript j denotes the site number 
and subscripts 7 and 8 correspond to F7 and Fg respec- 
tively. Using these states, we can derive the effective 
Hamiltonian H — X)™ '^"'^n'^"+i + + ^ncl^Cn, where 
c„/cjj denotes the renormalized annihilation/creation op- 
erator and the effective transfer integral T„ and site en- 
ergy Vn are given as: 
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FIG. 3: Localization length ^ of the effective single band 
model for Fg orbital with (a) <^ = 0°, (b) = 30°, (c) tp — 60°, 
and (d) (p — 90°, respectively. 
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can write down T, 
n,n+i is the angle between 



As for the Fg case, on the other hand, the effective 
site energy Vn explicitly depends on the local spin angle 
9n ■ If we have the spin configuration in the plane which 
is parallel to the zx plane, i.e., 9n =const., and set the 
pitch A(j) =const., Vn and T„ are constant. On the other 

^ 1 \ 

] — j hand, if we have the tilt of the spin rotation plane from 
the above plane to the other plane, 9n is no longer a 
— — ), constant and then Vn would generally be IC. |T„| also 
depends on both the angles of 5„ and Sn+i- Here we 
can conclude that the case where holes are in Fg, the 
as effective ID model would generally be IC. This explains 



why the upper part of the t2g density of states in Fig. 
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1 are localized more strongly, where the wavefunction is 
mainly from Fg components. 

Now we focus on the Fg case and numerically examine 
whether the localization of the wavefunction occurs in 
more details. We consider the helical spin configuration 
Sn = {S cos{qn) , S COS (fi sui{qn) , S sin ip sin(gn)), where 
q is helical wavenumber, and (p denotes the tilt angle of 
the spin rotating plane from a;2;-plane (See Fig. EJ. 

The numerical calculations are performed for systems 
of size 1009, a large prime number, with nearly incom- 
mensurate modulations q/2TT = j'/lOOQ, {j = 1, 2, 3, ...). 
The results are shown in Figs. 3, where the vertical and 
the horizontal axes represent the energy and the helical 
wavenumber, respectively. We take the unit where t = 1 
and a = 1. The tih angles are (p = 0°,3G°,60° and 90° 
for Figs. 3(a),(b),(c), and (d), respectively. The energy 
spectrum in Fig. 3(d) is almost same as the lowest band 
of Fg bands in Fig. 1. In Figs. 3, the localization length 
increases as the color changes from blue to green. The 
figures clearly display that there are domains of strong 
localization ^ ~ 1 when we have a finite tilt angle ip. On 
the other hand, for ip = 0° the transfer integrals r„ and 
the on-site potentials Vn are constants, and there is no 
localized states. Even in the most suitable case for local- 
ization, i.e.. Fig. 3(d), however, the helical wavenumber q 
should be approximately in the range of 27r/3 < q < An/ 3 
for the localized states. This is because the long period of 
the helical structure means the slowly varying and weak 
perturbations in the rotated frame, and hence does not 
cause the localization. 

Now we discuss the possible realization of the localized 
states in realistic systems. From the above results, three 
important conditions for the localization are (i) strong 
SOI, (ii) short helical period, and (iii) the direction of 
rotating spin plane. The SOI increases as the mass of the 
atom get heavier, and hence the present model becomes 
more relevant from the viewpoint (i). For 3ci-orbitals of 
transition metal oxides, the SOI is typically of the or- 
der of 20-30meV, which is an order of magnitude smaller 
than the transfer integral t. Therefore the localization 
length is expected to be rather large, and hence the dis- 
order efi^ect such as the impurity scattering mi ght hide 
the IC effect. Therefore, even though /3-Mn02 [ll| and 
SrFeOg (IlEl show interesting transport properties, it 
is unlikely that the localization found in this paper is 
relevant to these materials. 4d or 4/, 5/-orbitals, where 
SOI is larger than ^ 0.3eV, are more promising. Actu- 
ally there are many rare-earth metals sho wing the helical 
spin structure such as Tb, Dy, Ho, Er [2^l28l|. From the 
condition (ii), it is rather hard to find the short period 
helical spin structure. It is typically 4a-5a or even larger 
in 111 IH El. From this viewpoint, MnOa is an 

interesting case, but the localization is unlikely as dis- 
cussed above. As for the condition (iii), we need more 
study since only the cubic case has been considered. The 
directional dependence of the spin plane might be useful 



to control the localization/delocalization by an external 
magnetic field. 

Even though the conditions for localization discussed 
above are rather stringent, which explains why it has 
never been observed experimentally thus far, it will play 
a vital role in the quantum transport properties of the 
system once realized. One direct consequence is the large 
anisotropy of the resistivity between parallel and perpen- 
dicular to the helical axis, i.e., it should be much more 
resistive in the parallel direction. 

In conclusions, we have studied the localiza- 
tion/delocalization of the electronic states in helical mag- 
nets. We found the localized states under the condition 
of (i) strong spin-orbit interaction, (ii) short helical wave- 
length, and (iii) proper direction of the plane on which 
spins rotate. The strong dependence of the localization 
length ^ on the orbital is also found, which is explained 
by an effective model for a certain limiting case. 
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